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ABSTRACT 

We present 44 days of high duty cycle, ultra precise photometry of the 13th 
magnitude star Kepler-5 (KIC 8191672, T eff =6300 K, log 0=4.1), which exhibits 
periodic transits with a depth of 0.7%. Detailed modeling of the transit is con- 
sistent with a planetary companion with an orbital period of 3.548460 ±0.000032 
days and a radius of 1.43li°)!o52 -^J- Follow-up radial velocity measurements with 
the Keck HIRES spectrograph on 9 separate nights demonstrate that the planet 
is more than twice as massive as Jupiter with a mass of 2.114^° ) '° ) 5 ) 6 , Mj and a 
mean density of 0.894 ± 0.079 g/cm 3 . 
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Subject headings: planetary systems — stars: individual (Kepler-5, KIC 8191672, 
2MASS 19573768+4402061) — techniques: spectroscopic — Facilities: The Ke- 
pler Mission. 



1. INTRODUCTION 

The launch of Kepler offers a special opportunity to study the nature of transiting ex- 
trasolar planets through transit photometry. The Kepler mission is designed to determine 
the frequency of terrestrial-size planets. The Kepler mission also provides an excellent plat- 
form to study and understand the wide variety of extrasolar planets. From the spacecraft's 
Earth-trailing orbit heliocentric photometric measurements are devoid of artifacts usually 
associated with instruments in close proximity to the Earth, due to problems such as the 
day-night cycle and effects imposed by the Earth's atmosphere on ground based observations, 
and orbital effects for satellites in low earth orbit. 

We report the discovery and confirmation of Kepler-5b, a strongly irradiated transiting 
hot Jupiter. We describe the Kepler photometry and transit modeling, and the follow-up 
observations used to confirm that Kepler-5b is a planet, including an orbital solution using 
radial velocities obtained with HIRES on Keck 1. 



KEPLER PHOTOMETRY 



The Kepler photometer consists of a 0.95-m Schm idt telescope feed ing an array of 42 
CCDs with a broad (430 - 890 nm) spectral response (jKoch et al.ll2010l ). From its Earth- 
trailing heliocentric orbit and single field of view, Kepler can monitor more than 100,000 
stars nearly continuously for the lifetime of the mission. 

Photometry of Kepler-5b (KIC 8191672, a = 19 h 57 m 37 s .68, 5 = +44°02'06".2, Kep- 
mag=13.369) was obtained during 1 May - 15 June 2009, at an integration time of 29.426 
minutes for each observation. This run includes 9.7 days of photometric data taken while 
ground analysis was being pe rformed in prepa ration for start of science operations, and is 



referred to as the Q0 data set (IHaas et al. 



20101 ). Kelper science operations began on 12 May 



2009, which compose the second data set of 33.5 days of photometry referred to as the Ql 
data set. There is a gap with a length of 30 hours between Q0 and Ql to facilitate contact 
between the spacecraft and groundstations for downlink of the Q0 data stream. 



After pipeline data processin g and photometry extraction, as described in I Jenkins et al. 
(120101 ) and ICaldwell et al.l (120101 ) , the time series stream was detrended with a running 1-day 
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mean. All observations that occurred during a transit were rejected from the evaluation of 
the mean. The detrended photometry is presented in Figure [TJ 



2.1. Initial Transit Fits and False-Positive Detection 



From examination of QO light curves, Kepler-5 and a few dozen other candidates were 
initially flagged as Kepler Objects of Interest (KOI). To identify the best candidates for 
ground-based follow-up observations, a number of metrics are computed to help recognize 
and reject stellar binary systems that can mimic a planetary transit signal. A quick summary 
of these steps, as it relates to Kepl er-5 are listed below . For a full description of the steps 
taken to identify false positives see iBatalha et al.l (120101 ) . 



The transit lightcurve is modeled using the analytic express i ons o f iMandel fc Ago! 



(120021 ) using V-band non-linear limb darkening parameters IClaretl (120001 ) . The stellar ra- 
dius (.R*), effective temperature (T e ff) and surface gravity (logg) are fixed to values adopted 
from Kepler Input Catalog (KIC). The stellar mass (M*) is usually calculated from log g 
and R+. In the case of Kepler-5, a larger stellar radius of 1.9 R Q was required to be con- 
sistent with the transit duration. With M* and R* fixed to their initial values, a transit fit 
is then computed to determine the orbital inclination, planetary radius, and depth of the 
occultation (passing behind the star) assuming a circular orbit. T he best fit is found using 
a Levenberg-Marquardt minimization algorithm (IPress et al.lll992l ). 



These initial fits are used to determine whether the transit is representitive of a planetary 
event. In the case of a stellar binary where the surface brightness of the two components 
differ, the observed depth of the odd and even numbered transits will differ. There was no 
measurable difference in the transit depths to report. 

For systems that show only primary transits, a search is made for a weak occultation 
assuming that the orbit is circular, namely at phase 0.5. In the cases where an occultation 
is found with significance greater than 2cr for the depth, the dayside temperature of the 
planet can be estimated. From the depth of the occultation, the flux ratio of the planet 
and star (Fp/F+) over the instrumental bandpass can be obtained. The depth of the transit 
indicates the ratio of the planet and star radii. By assuming that the star and planet both 
behave as blackbodies and the flux ratio is bolometric, the dayside effective temperature can 
be estimated. For Kepler-5 we find Fp/F* = 3.6 x 10 -5 . This gives a planetary effective 
temperature of T e g = 1720 ± 214 K, where an error of 30% is assumed for the input stellar 
luminosity and radius. This estimate is a lower limit, as a significant fraction of the planetary 
flux is emitted at wavelengths longer than the red edge of the Kepler bandpass, but it is a 
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useful diagnostic in order to determine whether the depth of the occultation is consistent 
with a strongly irradiated planet. To make this comparison we can estimate the equilibrium 
temperature, 

T C(l = T^Rj2a)^[f(l-A B )f\ (1) 

for the companion, where R* and T* are the stellar radius and temperature, with the planet 
at distance a with a Bond albedo of As, and / is a p roxy for atmosphe ric thermal circulation. 
We assume A^ = 0.1 for highly irradiated planets (IRowe et al.ll2006l ) and / = 1 for efficient 
heat distribution to the night side. Assuming stellar irradiation is the primary energy source, 
these parameter choices give a rough estimate for the dayside temperature of the planet. 
Assuming a 30% error in the input stellar parameters and that the star and planet act as 
blackbodies we find T eq = 1810 ± 289 for Kepler-5b. The consistency of T e g and T eq to first 
order suggests that the occultation is consistent with a strongly irradiated planet. If the 
estimate of T e g- were found to be much larger than T eq , then the companion is likely to be 
self-luminous and is probably a star. 



2.2. Centroid Shifts 

The Kepler pipeline provides measurements of the centroids of stellar images with a 
precision of about 0.1 millipixel for the brightness of Kepler-5 (Batalha, 2010). Examination 
of the image centroids during transits revealed shifts with an amplitude of 0.2 millipixels. 
Such a motion can result if the photo-aperture for Kepler-5 includes light from another source. 
We examine and demonstrate with speckle and adaptive-optic (AO) images in Section 13.11 
that the centroid motion observed for Kepler-5 can be accounted for by the presence of 
nearby faint stars with constant brightness. 



FOLLOW-UP OBSERVATIONS 



After a KOI passes t he above tests , addit ional ground-based follow-up observations are 
obtained as described in iGautier et al.l (120101 ). These observations include high- resolution 
imaging to search for additional sources of flux within the Kepler photometric aperture that 
would dilute the depth of the transit, and reconnaissance spectroscopy to confirm and refine 
the KIC stellar classification and to search for evidence of stellar companions. 
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3.1. High-Resolution Imaging 

Ground-based visible-light speckle imaging from the WIYN Telescope and near-infrared 
AO imaging from the Mt. Palomar 5-m telescope show that Kepler-5 has two companions. 
The first companion is 0.9" away and 5.2 magnitudes fainter, as seen in the NIR-AO image, 
but is not visible in the speckle image. The other companion, KIC 8191680, is 7.3" away 
with a Kepler magnitude of 17.69. The difference in Kepler magnitude between Kepler-5 
and KIC 8191680 can be used to compute the expected centroid shift, depending on which 
object is assumed to be the source of the transit signal. 

If KIC 8191680 is a background eclipsing binary, then the centroids are expected to shift 
by +14.0 and —0.8 millipixels in the column and row directions of the detector, respectively. 
If KIC 8191680 shows no flux changes and Kepler-5 is indeed the source of the transit 
signal, then the expected centroid shift is —0.3 and 0.0 millipixels, which is consistent with 
the centroid shift reported in Section 12.21 The transit depth is estimated to be diluted by 
~ 2 ± 0.2%. We include the effects of dilution in our transitfits. 

The speckle image shows a single star within its 2" square field. This rules out very 
close background eclipsing binaries that might simulate the observed transits, except for 
stars lying closer than about ~ 0.1" from the target star. 



3.2. Reconnaissance Spectroscopy 

The Fibre-fed Echelle Spectrograph (FIES) on the 2. 5-m Nordic Optical Telescope 
(NOT) was used on 4 June 2009 to obtain a spectrum of Kepler-5 for the purposes of 
stellar classification. Stellar parameters of logg = 4.0 and T e g = 6500 K were estimated, 
which were used to refine the transit model fits described in §2.11 The spectrum shows no 
indication that the system includes an eclipsing binary or has a composite spectrum. 



3.3. HIRES Spectroscopy 

HIRES spectra were obtained from 3-6 June, 2-4 July, and on 6 October, 2009. Fig- 
ure [2] shows the radial velocities for Kepler-5 folded with the photometric period of the 
planet. An analysis of t he Keck/HIRES template spectrum by D. Fischer using SME 
JValenti fc Piskunov! fl~9~9~6l ) measured the stellar parameters as listed in Table [2j A fit to 
the velocities with the eccentricity fixed to zero measures a reflex velocity semi-amplitude 
of K = 228 ± 8 ms -1 with residuals of 15 ms" 1 . A bisector analysis of the HIRES velocities 
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does not show any significant variations correlated with the radial velocities. 



4. ANALYSIS 



The Kepler photometric and HIRES radial- velocity measurements can be simultaneously 
modeled. The data is fit for the center of transit time, period, impact parameter (6), the 
scaled planetary radius (R p /R*), the amplitude of the radial velocity (K), photometric and 
velocity zeropoints and £/R*. The last t erm, C/R* is rel ated to the transit duration (T^ = 
2((/R ir )~ 1 ) and the mean stellar density (jPal et al.l 120091 ) . We also account for non-circular 
orbits by modeling for e cos(o;), e sin(o;), where e is the orbital eccentricity and u is argument 
of periastron. 

The modeled stellar density is strongly dependent on the impact parameter, b, of the 
planetary orbit. It was discovered that the adopted limb-darkening parameters were placing 
unrealistic constraints on b. Specifically, the best \ 2 would occur at b = for a majority 
of the models of Kepler light curves. The model fits were repeated for various choices of 
limb darkening, including values generated specifically for the Kepler bandpass with Atlas 
9 models by A. Prsa. Detailed examination of the fits revealed that the limb darkening 
models were over-predicting the curvature of the base of the transit shape. To alleviate 
this problem, limb darkening coefficients wer e computed for the t hree known exoplanet 



systems in the Kepler fiel d of view, TrES-2 (IDaemgen et al.l l2009l ). HAT-P-7 (jPal et al. 



20081 : iGilliland et al.l 120101 ). and HAT-P-11 flBakos et al.ll2009f ). based on fits fixed to the 
published values of M*, R+, and i. We then linearly interpolate between these values for 
candidates with different temperatures. It is estimated this procedure results in a ~ 0.5° 
systematic error on the quoted value of i reported in Table [2J 

The error distribution for the stellar densities that fit the transit light curve are obtained 
from a Markov-Ch ain- Monte-Carlo (MCMC) analysis. We adopt the approach outlined in 
§4.3 of iFordl ( 120051 ). Convergence was tested by generating 10 Markov-chains each with 10 6 
steps and different initial conditions. We checked that each chain had statistically similar 
means and distributions. Given the stellar effective temperature from the Keck/HIRES tem- 
plate and the mean stellar density from transit photometry we employ the p* method. This 
method uses stellar evolution Yale-Yonsei tracks to calculate the range of stellar parame- 
ters consistent with the observations. A full description of the p* method is described in 
citetborlOb. 



The output of the p* method is a set of Markov Chains giving a consistent set of 10 5 
allowed M* and R± pairs. The distribution of stellar mass and radius allows us to define 
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a transition probability to compute a MCMC analysis to determine the most likely model 
values. We fit for stellar mass and radius, planetary mass and radius, center time of transit, 
orbital period and inclination, depth of the occulation and allow for non-circular orbits 
through the inclusion of e cos(o;) and e sin (a;). We generated 10 sets of Markov-Chains with 
10 6 elements and use all 10 sets to generate distributions for the model fits. 

The errors in Table [2] are centered around the mode of each model parameter which 
represents the value of highest relative probability. The mode is calculated by binning the 
bootstrapped parameters into 50 equally-spaced bins. We chose to use the mode to determine 
the most likely system parameters as some distributions, such as the stellar mass, are bi- 
modal. We report the bounds that enclose ±68% of the samples centred on the mode and 
report these values in Table [2J 



DISCUSSION 



Similar to Kepler-7 (ILatham et al.l 120101 ) . the Kepler-5 host star is not much hotter 
than the Sun, T e g = 6297 ± 60, but is much more massive and larger than the Sun, M* = 



1 Q7/I +0.040 



M© and R* = 1.793+^ R Q . With a mean stellar density of = 0.33 ± 0.02 
gcm~ 3 the star is on its way to becoming a shell-burning subgiant. Similar to Kepler- ^ 



(IBorucki et al.l l2010bl ) and Kepler-7 there are two distinct sets of model parameters that 
fit the observations equally well. There are two peaks in the mass distribution centered at 
M* = 1.21 and 1.38 M . While our analysis favours the larger mass, the lower mass choice 
can not be ruled out. 



The planet is found to have a mass of Mp 

I +0.041 



2.114±£°| Mj and a radius of R P 



1.431^052 -Rj, where Rj refers to the equatorial radius of Jupiter. This makes the planet 
more massive and larger than Jupiter. The mean planetary density is 0.89 ± 0.08 gem -3 , 
which is not u ncommon amo n gst th e known population of transiting Jupiter-sized planets. 
Figure 3 from ILatham et al.l (120101 ) compares the mass and radius of Kepler-5b to other 
known transiting extrasolar planets. 

The photometry points to a weak (2a) detection of an occultation. As correlations in the 
timeseries may lead to a false detection of an occultation we computed a Fourier Transform 
and autocorrelation of the dataset which did not show any significant issues. There is no 
excess of power in the Fourier Transform from 1 to 12 hours has a 3 a detection limit of ~15 
ppm. 



Most of the observed light from the planet in the Kepler bandpass is due to thermal 
emission, but we can use the depth of the occultation to derive an upper limit on the albedo 
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of the planet. The geometric albedo, 

F a 2 

A 9 = (2) 

is defined as the ratio of the planet's luminosity at full phase to the luminosity from a 
Lambertian disk, where R p is the planetary radius. This allows us to place boundaries on 
the dayside temperature of the planet as estimated by the equilibrium temperature given 
by Equation [TJ Figure [3] shows the boundaries for Kepler-4b, 5b, 6b and 7b where we 
have calculated T eq for / =1,2 and A-g from zero to the 1-cr upper limit. The Bond albedo 
is estimated from the geometric albedo by assuming a Lambertian scatter as described in 



Rowe et al.l (120061 ). Due to the relatively short periods of the first four planets discovered by 
Kepler, it not surprising that all of them have temperatures greater than 1500 K. The upper 
limits show that Kepler-5b, 6b and 7b are less reflective than Jupiter and likely quiet dark, 
consistent with predictions and observations of other hot extrasolar planets. Additional 
photometry of these planets will provide special opportunities to study planets in extreme 
conditions. 



Funding for this Discovery mission is provided by NASA's Science Mission Directorate. 
We would like to thank D. Fischer, A. Prsa and everyone that has contributed to the Kepler 
mission. 
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Fig. 1. — The phased light curve of Kepler-5 containing 12 transits observed by Kepler 
between 1 May and 15 June 2009. The upper panel shows the full 44-day time series after 
detrending. The bottom panel shows the light curve folded with the orbital period. The 
lower curve shows the primary eclipse, with the fitted transit model overplotted in red and 
corresponding scale to the left. The upper curve covers the expected time of occultation, with 
the fitted model overplotted in green and corresponding scale found to the right, expanded 
by a factor of 5 relative to the scale for the transit data. 
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Fig. 2. — The phased radial- velocity curve for Kepler-5 consisting of 8 epochs observed 
using the Keck/HIRES spectrometer spanning 126 days. The overplotted fit assumes a 
circular orbit, phased to match the transit photometry. A tabulation of the radial-velocity 
measurements used for this Figure may be found in the online materials. 
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Fig. 3. — The dayside temperatures and 1 a upper bounds on the Bond albedo for Kepler-4b 
(hatched), Kepler-5b (solid), Kepler-6b (outside) and Kepler-lb (cross hatched) based on 
the 1-a upper limits on the the depth of the occultation. 
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Table 2. System Parameters for Kepler-5 



Parameter 

Orbital period P (d) 
Midtransit time E (HJD) 
Scaled scmimajor axis a/Ft* 
Scaled planet radius Rp/R± 
Impact parameter b = a cos i/i?* 
Orbital inclination i (deg) 
Orbital semi-amplitude K (ms -1 ) 
Orbital eccentricity e 
Ccntcr-of-mass velocity 7 (ms -1 ) 

Observed stellar parameters 
Effective temperature T e g (K) 
Spectroscopic gravity logg (cgs) 
Metallicity [Fe/H] 
Projected rotation vsini (kms -1 ) 
Mean radial velocity (kms -1 ) 

Derived stellar parameters 
Mass M*(M Q ) 
Radius R*(Rq) 
Surface gravity logg* (cgs) 
Luminosity (Lq) 
Age (Gyr) 

Planetary parameters 
Mass M P (M,j) 
Radius Rp (Rj, equatorial) 
Density pp (gem -3 ) 
Surface gravity loggp (cgs) 
Orbital scmimajor axis a (AU) 
Equilibrium temperature T eq (K) 



Value Notes 



3.548460 ± 0.000032 A 

2454955.90122 ± 0.00021 A 

6.06 ±0.14 A 
0.081951°" A 

0.3931HS A 

86?3 ± 0.5 A 

227.5 ±2.8 A,B 

< 0.024 A,B,G 

A,B 

6297 ±60 C 

3.96 ±0.10 C 

±0.04 ± 0.06 C 

4.8 ±1.0 C 

-46.7 ±4.1 B 

1.374t°;g*§ C,D 

1.793l°;gg CD 

4.07 ±0.02 C,D 
4.67l°;g C,D 
3.0 ±0.6 C,D 

2.114±°;ggg A,B,C,D 

L431 io:o52 A,B,C,D 

0.894 ±0.079 A,B,C,D 

3.41 ±0.03 A,B,C,D 

0.05064 ± 0.00070 E 

1868 ± 284 F 



Note. - 
A: Based on the photometry. 
B: Based on the radial velocities. 
C: Based on a MOOG analysis of the FIES spectra. 
D: Based on the Yalc-Yonsei stellar evolution tracks. 
E: Based on Newton's version of Kepler's Third Law and total mass. 
F: Assumes Bond albedo = 0.1 and complete redistribution. 
G: 1 sigma upper limit 



